7* ^>T7T7r^  7?  '.^ '.v  ;v  ■■v ».  ■:  ■  t  _  * ,  » ■.  i  ■  ^  ■  1; 


\  v.  -...%*:<  .  .v  **.■■/■/ 


is 


is 

fe 

Sv 

tv 


SP* 

i 


Cu 


<1 

§ 

No 


<3 

\ 


INVESTIGATION  OF  METALLIC  IMPURITIES 
INTRODUCED  INTO  SI02 

Mid-Term  Technical  Report 
For  Period  9/1/82  through  9/31  /83~~j 
Contract  No.  MDA903-82-k-0412 
DARPA  Order  No.  AO  4541 

Effective  Date  of  Contract:  06/6/82 
Expiration  Date  of  Contract:  05/31/84 

Principle  Investigators: 

R.  M.  Swanson  and  T.  W.  Sigmon 
(415)  497-0426 

/  ' 

|S*  • 

Stanford  University  V' 


Sponsored  by 


JEC  1  9  1933 

r' 

A 


Defense  Advanced  Research  Projects  Agency  (DOD)  DARPA  Order  No.  AO 
4541  under  contract  No.  MDA  903-82-K-0412  issued  by  Department  of  Army, 


Defense  Supply  Service- Washington,  Washington,  D.  C.  20310. 

Disclaimer 

The  views  and  conclusions  contained  in  this  document  arc  those  of  the  authors 
and  should  not  be  interpreted  as  representing  the  official  policies,  cither  expressed 
or  implied,  of  the  Defense  Advanced  Research  Projects  Agency  or  the  U.  S. 
Government. 

approved  for  public  release 
distribution  UNLIMITED 


ft 


88  10  2'o  9  0 


SUMMARY 


---^This  report  contains  the  theoretical  development  of  a  model  that  describes  the 
diffusion  of  interstitial  diffusors  in  SiO <£•  Experimental  data  supporting  this  model 
and  the  implications  and  conclusions  predicted  by  this  model  are  presented. 

In  summary  this  model  concludes  that  Pd,  Au,  Mo,  Ta,  W,  Pt,  Ti,  and  A1  are 
potential  candidates  for  use  in  VLSI  gate  metalization  schemes,  where  as  Ag,  Cu, 
Na,  Ni,  Mn,  Fe,  Mg  and  Ga  are  unsuitable.  (r~ 
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Chapter  I 
MODELING 


A  complete  derivation  of  the  diffusion  equation  is  found  in  the  Appendix  along 
with  analytical  solutions  for  some  interesting  boundary  conditions.  In  this  chapter 
a  new  diffusion  and  solid  solubility  model  is  presented  that  allows  calculation  of 
the  diffusion  and  solid  solubility  activation  energy.  Comparison  of  this  model  with 
experimental  data  is  done  in  chapter  II. 

A.  Interstitial  Diffusion 

The  diffusion  of  metals  in  silicon  dioxide  is  assumed  to  be  interstitial  in  nature 
[1.1,  1.2]  .  The  diffusivity  or  diffusion  coefficient  for  this  type  of  diffusion  is  given 
by  the  following  expression: 

D  =  Dae TH*  (1.1) 

where  Ea  is  the  activation  energy  for  diffusion,  k  is  the  Boltzmann  constant,  T  is 
the  temperature  and  Da  is  the  diffusion  constant. 

To  evaluate  the  diffusion  coefficient  at  a  given  temperature  the  values  of  D0 
and  Ea  must  be  known.  The  following  two  sections  will  allow  us  to  estimate  their 
values  in  terms  of  fundamental  quantities. 

1.  Diffusion  Constant,  D0 

Using  the  Zener  theory  [  1.1  ]  to  predict  D0  allows  us  to  write 

D0  =  pb2vr)fe~£  (1.2) 


W 


-  1  - 


where  p  is  the  probability  of  an  interstitial  jump,  b  is  the  jump  distance,  rj  is  the 
number  of  equivalent  jump  paths,  u  is  the  lattice  vibration  frequency,  /  is  the 
correlation  factor,  k  is  the  Boltzmann  constant  and  Sa  is  the  entropy. 

From  this  equation  we  see  that  D0  depends  upon  the  crystal  structure  of 
the  lattice  through  which  the  dilfusor  is  traveling.  In  all  cases  Sa  will  be  a  positive 
number  [  1.1  ]  so  that  D0  will  have  a  lower  limit  of 


D0  >  pb2urjf 

To  first  order 


(1.3) 


So  that 


prjf^l 


(1.4) 


D0  b2u  .025 cm2 /sec  (1.5) 

is  the  lower  limit  for  D0.  This  calculation  used  a  lattice  spacing  for  SiO-2  of  5  A 
and  a  lattice  vibration  frequency  of  1013hz/sec  [1.3,  1.4]  We  will  use  this  number 
as  an  estimate  of  Da  in  our  future  predictions  of  D. 


2.  Activation  Energy  of  Diffusion,  Ea 

To  calculate  the  activation  energy  of  metal  diffusion  in  silicon  dioxide  we 
use  a  classical  approach  introduced  by  Anderson  and  Stuart  in  1954  [1.5].  In  this 
model  they  assume  that  the  activation  energy  is  the  sum  of  two  terms;  the  first 
is  the  strain  energy  between  the  lattice  structure  and  the  diffusing  atom,  and  the 
second  is  the  electrostatic  energy  between  a  charged  diffusing  species  and  the  silicon 
dioxide  lattice.  We  consider  each  in  the  following  sections. 

-  2  - 


a.  Strain  Energy 


We  first  assume  that  the  strain  energy  generated  in  a  silicon  dioxide 
network  due  to  the  enlargement  of  a  spherical  doorway,  of  radius  rd,  to  accom¬ 
modate  an  ion  of  radius  r  is  approximately  the  elastic  energy  required  to  dilate  a 
spherical  cavity  from  radius  rd  to  r.  This  strain  energy  (Es)  is  given  by 

Ea  =  47T Grd(r  -  rd)2  (1.6) 

where  G  is  the  shear  modulus  for  silicon  dioxide,  G  =  3.12  X  1011  dynes /cm2  [1.4], 
rd  is  the  radius  of  the  window  through  which  the  diffusing  atom  must  pass,  rd  = 
0.6  A  [1.5],  and  r  is  the  radius  of  the  diffusing  atom.  The  above  expression  is 
reduced  by  a  factor  of  two  from  the  exact  mathematical  solution  to  account  for  the 
very  loose  structure  of  silicon  dioxide. 


b.  Electrostatic  Energy, 


We  will  approximate  the  electrostatic  energy  change  by 

Eb  _  Pzz0e2 


(1.7) 


7  (r  +  r0) 

after  Anderson  and  Struart  [1.5].  where  /3  is  the  finite  displacement  factor  in 
angstroms,  given  by 

2.1  -r 


/?  = 


3.5 


(1.8) 


and  r  is  the  radius  of  the  diffusing  species.  The  geometric  factors  arc  determined 
by  assuming  an  interstitial  diffusion  jump  distance  of  7  A  and  an  atomic  radius 
for  oxygen  (r0)  of  1.4  A.  The  charge  states  of  the  diffusing  species  and  the  oxygen 
in  the  silicon  dioxide  lattice  are  z  and  z0  respectively,  for  oxygen  z0  is  2.  7  is 
approximately  the  dielectric  constant  of  silicon  dioxide,  7  =  3.9. 


ACTIVATION  ENERGY  (eV/PARTICLE 


RADII  (A) 


Figure  1.1.  Interstitial  Diffusion  Activation  Energy.  Interstitial  diffusion  activa¬ 
tion  energy  calculated  by  the  sum  of  the  strain  energy  and  the  change 
in  the  electrostatic  energy. 
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c.  Activation  Energy,  Ea 

Combining  the  strain  energy,  E3,  and  the  electrostatic  interaction 
energy,  E 5,  yields  the  total  interstitial  diffusion  activation  energy.  The  material 
parameter  table  found  in  Appendix  A  contains  the  known  atomic  radii  for  various 
neutral  and  ionized  metal  atoms.  Using  these  numbers  and  those  presented  allows 
the  generation  of  Fig.  [  1.1  ].  Figure  [  1.1  ]  allows  us  to  predict  the  activation  energy 
for  interstitial  diffusion  in  S1O2  for  atoms  with  known  atomic  radii. 

3.  Diffusion  Constant  for  Interstitial  Diffusion 

Having  values  for  D0  and  Ea  allows  prediction,  via  Eqn.  (1.1)  of  the 
interstitial  diffusion  coefficient  for  elements  that  we  believe  diffuse  interstitially. 
Figures  [1.2  —  1.4]  present  the  calculated  Diffusivity  for  various  metals  and  their 
ions  in  SiO^.  The  data  used  in  the  calculation  of  the  figures  arc  found  in  Appendix 
A.  Values  appearing  in  the  tables  in  brackets  represent  estimated  values,  for  ele¬ 
ments  which  no  published  data  could  be  found. 

In  the  next  section  we  derive  an  equation  for  determin:ng  the  solid 
solubility  for  these  metals  in  silicon  dioxide. 
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1.2.  Predicted  Neutral  Diffusion  Coefficients.  Calculated  Diffusivity  for  the 
indicated  metals  in  silicon  dioxide  as  a  function  of  temperature. 
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Figure  1.4.  Predicted  Ionized  Diffusivity.  Calculated  Diffusivity  for  doubly  and 
triply  ionized  atoms. 


B.  Solid  Solubility 


The  solid  solubility  of  metals  in  silicon  dioxide  is  derived  from  solutions  of  the 
thermodynamic  equations  for  each  possible  diffusing  species.  See  Appendix  B  for 
derivation  of  the  law  of  mass  action  and  chemical  potential. 

The  important  assumption  that  makes  this  derivation  justifiable  is  that  the 
ideal  gas  law  holds  for  the  loose  structure  of  silicon  dioxide.  This  allows  us  to  calcu¬ 
late  the  concentration  of  metal  in  silicon  dioxide  as  a  function  of  the  concentration 
of  metal  vapor.  In  the  following  sections  we  will  derive  the  concentration  of  neutral, 
singly,  doubly,  and  triply  ionized  metal  atoms  in  silicon  dioxide  or  any  other  loosely 
bound  lattice  utilizing  the  above  assumption. 


1.  Neutral  Diffusors 

Using  the  law  of  mass  action  we  write  the  concentration  of  neutral  metal 
atoms  [Mox\  in  silicon  dioxide  as 

\Mox\  —  [Mgas]e  (1*9) 

where  E{nt  is  the  total  internal  energy  and  is  the  sum  of  the  internal  energy  of  the 
silicon  dioxide,  Emox,  plus  that  of  the  gas,  E\fg at.  The  internal  energy,  in  this  case, 
is  equal  to  zero  since  we  are  assuming  that  the  ideal  gas  law  holds,  i.e.  Emox  — 
EMgna.  k  is  the  Boltzmann  constant  and  T  is  the  temperature.  The  concentration 
of  metal  that  sublimes  from  the  solid  in  accordance  with  the  ideal  gas  law, 

,  can  be  expressed  as 


f  _ 1  = 


nTc-i\  _A 

- /r 


lnlO+(DTX  10_:'/nl0) 


(\  in'* 


In  the  above  equation  A,  B,  C,  and  D  are  constants  that  describe  the  vapor  pressures 
for  each  metal.  [1.6]  See  material  parameter  table  in  Appendix  A.  N  is  Avogadro’s 
number,  6.023  X  1023 atoms / mole.  R  is  the  gas  constant,  62396.0  torr  cm3  /  K  mole 
and  T  is  the  temperature  in  degrees  Kelvin. 

Thus  for  the  neutral  case,  the  amount  of  metal  entering  the  silicon  dioxide 
from  the  overlying  metal  film  is  that  predicted  by  the  ideal  gas  law.  Substituting 
the  above  constants  into  Eqn.  (L10)  and  rearranging  results  in 


—(Eint+AklnlQ—(kDT“  X  1 

[Mox]  =  9.653  X  10(fl+18)e— - & 


[kOT2  XlO~3lnlO)-kT(C-l)lnT) 


(1.11) 


which  is  plotted  in  Fig.  [  1.5  ]  for  various  metals. 


2.  Singly  Ionized  Diffusors 

To  find  the  solid  solubility  for  singly  ionized  diffusors  the  following  reac¬ 
tion  equations  must  be  solved: 


Mao[  Mgag  AHgub 


(1.12) 


M, 


gas 


A/'+  -4-  c 

lY1  gas  '  *gas 


A  E%on 


(1.13) 


M  ga3  M0%  &Eimage  d"  ^^binding  d”  ^^other  (1-14) 


e<jas 


ox 


^Eimage  d  X  +  A 2? other 


(1.15) 


cm 


TEMPERATURE  (°K) 


Figure  1.5.  Interstitial  Diffusion  of  Neutral  Atoms.  Solid  solubility  calculations  for 
the  interstitial  diffusion  of  neutral  atoms.  Pd,Mo,Fe,Au,Cu,Ta,VV,Ni 
and  Pt  have  solid  solubilities  less  than  lOI0cm~3. 


In  the  above  equations  the  heat  of  formation  for  each  reaction  is  shown  to  the  right 
of  each  equation.  A H3X1f,  is  the  heat  of  sublimation  and  is  known  for  all  elements, 
[1.3, 1.4].  A E{on  is  the  ionization  potential  and  is  also  known  for  most  elements, 
[1.3, 1.4].  A  ^binding  *s  found  as  described  in  section  26.  x  is  the  electron  affinity. 
A Eirnagc  *s  the  image  charge  for  both  the  positive  ion  and  the  negative  electron. 
AEother 's  any  other  heat  of  formation,  such  as  a  catalyst  to  the  diffusion  process. 

Combining  all  the  unknown  heats  of  formations  into  one  term, 
^remainder  >  gives 

remainder  ~  2AE image  4~  2AE0ther  4"  X  (1*16) 

Combining  the  above  reaction  equations  results  in 

M30l  M ox  +  eox  A Hgub  +  AEion  A E binding  A Eremainder  (1-17) 

for  the  metal  solid  solubility  for  the  singly  ionized  diffusor. 

We  can  deal  with  this  equation  in  a  manner  similar  to  the  neutral  case, 
which  gives  the  following  expression  for  the  solid  solubility  of  singly  ionized  atoms 
in  SiO-2  (See  Fig.  [1.6]). 


[M+ ]  =  1.527  X  10(£+17)e 


-(Eint  +  kAlniO-(kT2  DX  U)-3/nlO)- ££(2C+l)inT) 


1FT 


(1.18) 


3.  Doubly  Ionized  Diffusors 

For  doubly  ionized  diffusors  we  can  derive  the  following  reaction  equation 
30l  M ox  +  2e  A IIsut)  +  AEion  4-  2A -^remainder  (1.19). 

In  this  derivation  we  have  assumed  that  the  unknown  remainder  heat  of  formation 
is  twice  that  of  the  singly  ionized  case.  Applying  the  law  of  mass  action  in  a  manner 
similar  to  the  previous  two  cases  results  in 


i 


* 


[Mjos 


=  e 


— Xflos-Xw++-2xe- 


lMi+][=-]2 

where  x  for  each  species  is  defined  in  Appendix  D.  We  also  know  that 


(1.20) 


[M++]  =  2[«T] 


(1.21). 


Combining  these  equations  and  reducing  we  obtain 

f  2irmekT\  zEmt 


[m++] = 4  [Mgaa^(^^y-^ 

Where  E{nt  is  the  internal  energy  of  each  species 


(1.22). 


A Em  =  2AE,-  -  4£,„  -  A E„++ 


gas 


'Ml 


(1.23) 


or 


&int  —  AEion  &Ebinding  2A  Eremainder 

The  solid  solubility  for  the  doubly  ionized  case  is  then 

•  .  ...  -(Eint+kAlnl0-(kT*Dxl0~3lnlO)-kT[c+2)lnT) 

[M++]  =  6.084  X  10(lb+3)e —  - - sir - 


(1.24) 


(1.25) 


and  is  plotted  in  Fig.  [1.7 


4.  Triply  Ionized  Diffusors 


The  solid  solubility  of  the  triply  ionized  diffusors,  [MjI++],  is  found  using 
an  approach  similar  to  the  other  cases  and  results  in. 


r > i  .Lf 2vmekT\*  -Hint 
[M ox  }  —  J  e 


(1.26). 


Assuming  that  the  unknown  heat  of  formation  is  three  times  that  of  the  singly 
ionized  case.  We  obtain 


Afe'int  A^ion  '  &E binding  ^-'remainder 


(1.27) 
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Figure  1.8.  Triply  Ionized  Metal  Solid  Solubility  in  Si'02 •  Calculation  of  the  solid 
solubility  of  triply  ionized  metal  atoms.  Na  has  a  value  -<  1010em-3. 
No  other  atoms  were  calculated  due  to  lack  of  published  information. 


In  general  the  solid  solubility  takes  on  the  form  of 


^ effective 

C3ol  =  CQe  CT*  (1.30) 

for  all  cases.  A  negative  value  for  Eej  jective  implies  that  the  atom  in  question  would 
rather  be  in  the  oxide  than  out.  Values  for  Eej fecnve  for  various  elements  are  found 
in  Appendix  A. 


5.  Calibration  of  Solid  Solubility  Calculations 

The  figures  shown  in  this  chapter  were  based  on  a  value  for  AETemai-n(ier 
calculated  from  the  measured  solid  solubility  of  silver  in  silicon  dioxide.  This 
calculation  assumes  that  all  of  the  silver  in  the  oxide  was  put  there  while  in  a  singly 
ionized  state.  In  light  of  the  data  presented  earlier  which  showed  no  detectable 
diffusion  for  unbiased  samples  this  is  quite  reasonable. 

Using  an  upper  bound  ,  as  observed  with  RBS,  of  [Agox\  =  3  X  1019cm~3, 
we  find  AErematnder  =  8.87eV.  Using  this  value  in  the  following  equation 

Einl  +  kAlnlO  -  (kT2D  X  10~3inl0)  -  *£(2 C  +  1  )lnT 
A£«//«k„«  = - 5 - 5 - - — — - - - (1.31) 

we  obtain  a  value  for  AEe^ ]tcnve  of  approximately  0.5eV  for  singly  ionized  silver 
in  Si02- 

For  the  doubly  and  triply  ionized  cases  AEremain£/er  was  taken  to  be 
9eV  since  we  can’t  be  sure  of  the  exact  values  without  more  accurate  data  for  the 
ionized  radii,  image  charge,  etc.  However,  these  values  allow  us  to  see  trends  of 
the  solid  solubilities  in  silicon  dioxide  and  should  be  close  to  the  actual  values.  In 
the  following  chapter  we  compare  this  model  with  actual  experiments  to  test  its 
validity. 
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Chapter  II 

EXPERIMENTAL  ANALYSIS 


In  this  chapter  a  simplified  expression  for  the  diffusion  equation  which 
allows  extraction  of  the  diffusion  coefficient  from  experimental  data  is  developed. 
The  extracted  diffusion  coefficients  and  the  observed  activation  energies  are  then 
compared  with  the  values  predicted  by  the  theory  developed  in  chapter  I. 

A.  Solution  of  Diffusion  Equation 

A  complete  solution  to  the  diffusion  equation  is  given  in  Appendix  C  and  Ref. 
[2.1  ].  In  this  section  we  present  a  solution  to  the  diffusion  equation  that  holds  when 
the  electric  field,  i ,  dominates  the  diffusion  process  and  still  meets  the  following 
critria, 


£  <  —  ^TX  10 22V/cm  (2.1) 

qb 

Then  the  flux  can  be  written  as  (see  Appendix  C) 

J  =  qnCaol£  -  D6-f  (2.2). 

ox 

In  this  equation  q  is  the  charge  of  an  electron,  b  is  the  distance  between  potential 
barriers,  about  5  A.  n  is  the  mobility,  /j.  ==  D,  C  is  the  concentration  and  x  is 
the  distance  into  the  silicon  dioxide. 

When  the  electric  field  term  is  large  compared  to  the  concentration  gradient 
Eqn.  (2.2)  can  be  further  reduced  to 

J  =  IfDC'oi'  (2-3)- 
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In  the  following  this  equation  is  used  to  determine  the  diffusion 
the  other  parameters  are  known. 


1.  Silver  Diffusion 

Using  measured  data  on  silver  diffusion  in  silicon  dioxide  at  300°C  we 
derive  a  flux  for  silver  through  silicon  dioxide.  This  is  done  by  dividing  the  amount 
of  silver  observed  to  pile  up  at  the  silicon/silicon  dioxide  interface  during  a  BTS 
by  the  anneal  time.  In  Fig.  [2.1]  we  show  how  the  amount  of  silver  observed 
at  the  interface  in  atoms  /cm2  increases  with  time.  The  slope  of  this  plot  is  the 
flux  of  silver  through  the  silicon  dioxide.  The  solid  solubility  of  silver  in  silicon 
dioxide  at  300°  C  was  observed  with  RBS  to  be  3  X  1018cm-3  as  an  upper  bound. 
Combining  this  number  along  with  an  electric  field  of  5.2  X  104  V /cm  and  a  flux  of 
1.43  X  1012cm-2/sec,  see  Fig.  [2.1]  we  obtain  by  rearranging  Eqn.  (2.3) 


Da_q  =  4.5  X  10  l3cm2/sec 


at  300°  (7. 

Using  measured  silver  diffusion  data  in  silicon  dioxide  at  various  tempera¬ 
tures  Fig.  [2.2]  can  be  drawn.  The  slope  of  this  plot  is  the  activation  energy,  E  —  a, 
for  silver  diffusion  in  silicon  dioxide,  Ea  is  approximately  1.24eV. 
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Figure  2.1.  Silver  Pile  Up  at  the  SiC^/S*  Interface  versus  Time. 
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a.  Model  versus  Data 


The  model  developed  in  Chapter  I  expresses  the  diffusion  coefficient 
by  Eqn.  (1.1)  and  the  solid  solubility  by  Eqn.  (1.30).  Combining  these  two 
equations  with  Eqn.  (2.3)  we  write  the  flux  of  atoms  through  the  Si02  as 

Q  “  f  f  U6ton  +  ^  *olidsolubility) 

J  =  kT£DoC°e - *T -  (24) 

We  see  from  the  above  equation  that  the  observed  activation  energy  is  equal  to  the 
sum  of  the  activation  energies  defined  in  Chapter  I.  From  Chapter  I  the  predicted 
activation  energy  for  diffusion  is  1.3eV  and  the  predicted  solid  solubility  activation 
energy  is  0.5eV.  The  sum  of  1.8eV  is  to  be  compared  with  that  observed  of  1.24eV. 

Evaluating  Eqn.  ( 1.1 )  at  300°  C  for  silver  results  in 


D^'dicted  _  o  025e-tr  =  9.2  X  10~14cm2/see  (2.5) 

which  is  to  be  compared  with  our  observed  value  at  300°C  of  D^g  —  4.5  x 
10-13cm2/sec.  We  will  reserve  further  discussion  on  these  comparisons  for  the 
discussion  section. 


2.  Copper  Diffusion 

To  calculate  the  diffusion  coefficient  for  copper  we  use  the  model  developed 
in  Chapter  I  to  estimate  the  solid  solubility  of  singly  ionized  copper  in  silicon 
dioxide  at  450oC  to  be  1  X  1018cm~3.  This  assumption  is  necessary  since  the 
solid  solubility  of  copper  in  silicon  dioxide  is  below  the  detection  limit  of  the  RBS. 
This  number  should  be  a  high  estimate  and  therefore  will  give  us  the  worst  case 
answer.  This  approach  also  assumes  that  all  the  copper  in  the  SiO •>  is  singly  ionized. 
Using  the  observed  data  for  copper  we  can  make  plots  similar  to  Figs.  [2. 1, 2. 2] 
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dioxide,  sec  Figs.  [2. 3, 2. 4]. 


We  attribute  the  saturating  of  the  data  to  the  oxidation  of  the  copper 
during  BTS.  After  the  copper  has  oxidized  to  the  Cu/Si02  interface  we  lose  both 
adhesion  and  electrical  contact  resulting  in  loss  of  the  field  aided  diffusion.  Before 
oxidation  of  the  copper  dominates,  we  find  the  observed  activation  energy  to  be 
1.8eV  and  the  flux  of  copper  atoms  through  the  silicon  dioxide  at  450°C  to  be 
approximately  1.0  X  1013cm-2/sec.  Equation  [2.3]  can  be  used  to  calculate  the 
diffusivity  of  copper  in  silicon  dioxide  at  450°  C  with  an  electric  field  of  i  —  5.2  X 
104  V /cm,  resulting  in 


Dcu(450°C)  =  7^— F  =  1.2  X  I0_11cm2/sec 

k sol £ 


(2.6). 


The  result  of  Eqn.  ( 2.6 )  is  in  reasonable  agreement  with  the  value 
of  copper  diffusion  at  450°C  predicted  by  Eqn.  (1.1) 

Dpredicted  _  Q  025e^r  =  1.2  X  10"1Ocm2/sec  (2.7) 

Summing  the  diffusion  activation  energy,  1.2eV,  and  the  solid 
solubility  activation  energy,  0.646^,  results  in  the  total  activation  energy,  1.84eV, 
which  compares  nicely  with  the  observed  value  of  l.8eV. 


3.  Palladium,  Gold  and  Titanium  Diffusion 


The  diffusion  of  gold  and  palladium  in  silicon  dioxide  was  undetectable 
using  both  C-V  and  UBS  techniques.  This  was  expected  in  light  of  our  model  which 


prcdictes  a  very  low  solid  solubility  in  Si02  for  these  metals.  Thus  both  palladium 
and  gold  are  potential  candidates  for  low  temperature  processing. 

Titanium  was  found  to  react  with  silicon  dioxide  forming  a  Titanium 
oxide  compound  and  is  therefore  not  expected  to  be  an  interstitial  diffusor.  This 
agrees  with  what  was  observed  experimentally.  Also  since  Titanium  is  quite  reactive 
the  model  presented  earlier  is  not  expected  to  apply. 

B.  Discussion 

A  comparison  between  the  experimental  results  and  those  predicted  by  our 
model  is  quite  good.  Excellant  agreement  is  found  for  copper.  A  slight  discrepancy 
was  found  in  the  case  of  silver  in  that  the  theory  predicts  a  lower  diffusivity.  In 
light  of  experiments  that  revealed  the  catalytic  behavior  of  the  atmosphere  in  which 
the  experiment  takes  place,  we  propose  the  following  explanations: 

(1) The  strain,  component  of  our  model  could  easily  be  modified  by  the  presences 
of  a  molecule  or  atom  which  changes  the  elfective  window,  (r^)  through  which 
diffusion  occurs.  The  presence  of  such  molecules  or  atoms  could  also  have  an  effect 
on  the  electrostatic  energy.  If  they  were  charged  they  would  be  able  to  shield  or 
enhance  this  effect  and  change  the  activation  energy  from  that  predicted  by  our 
model. 

(2) These  molecules  or  atoms  could  also  modify  our  description  of  the  solid 
solubility  by  adding  another  heat  of  formation  into  the  equations.  If  we  consider 
the  silicon  dioxide  as  having  a  band  gap  then  the  addition  of  a  molecule  or  atom 
into  Si()2  changes  the  Fermi  energy.  This  change  in  the  Fermi  energy  may  lower 
the  ionization  energy  of  the  metal  atoms  and  therefore  change  the  solid  solubility. 

This  comparison  of  our  model  with  experimental  data  illustrates  the  usefulness 
of  this  approach.  More  work  is  needed  before  this  model  is  completely  quantitative. 
In  its  present  state,  however,  it  can  be  used  to  predict  which  metals  arc  likely 
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Chapter  III 

IMPLICATIONS  AND  CONCLUSIONS 


In  this  chapter  recommendations  for  metals  which  maybe  used  successfully  for 
VLSI  applications  are  presented.  Reasons  for  these  conclusions  are  also  presented. 
Futher  results  of  experiments,  not  yet  documented,  are  presented  in  support  of  these 
predictions. 


A.  Implications 

To  determine  which  metals  are  viable  candidates  for  VLSI  we  need  to  define 
a  criteria  to  compare  each  metal  against.  In  the  sections  that  follow  we  use  the 
number  of  metal  atoms  found  at  the  SiOz/Si  interface  or  an  equivalent  as  our 
means  of  comparison.  Further,  we  define  an  integrated  area  of  10locm~2  or  greater 
as  the  point  of  failure  for  a  metal  at  a  given  temperature. 

We  present  the  neutral,  singly,  doubly  and  triply  ionized  cases  for  most  of  the 
atoms  listed  in  Appendix  A.  For  some  atoms  it  is  obvious  that  the  present  approach 
doesn’t  apply.  Discussion  of  these  discrepancies  will  occur  later. 


1.  Neutral  Case 


For  the  neutral  case  we  solve  Eqn.  (2.2)  with  <f  =  0  and  the  following 
boundary  conditions: 


C(oo,  t)  =  0 


C(x,  0)  =  0. 

Where  C  is  the  concentration  and  Ca  is  the  solid  solubility. 


The  solution  is 


C(x,  t)  =  Ca  er  fc 


2{Dt)* 


where  Ca  is  given  by  Eqn.  (1.30),  D,  the  diffusion  coefficient,  is  given  by  Eqn. 
( 1.1 )  and  t  is  the  time  in  seconds. 

To  get  from  this  equation  to  our  definition  for  failure  we  must  integrate 
from  x  to  infinity.  Here  x  is  the  location  of  a  fictitious  silicon/silicon  dioxide 
interface.  This  gives  us 


CNt  — 


i 


eTUT  —  x  erfc 


2  (Dtp 


The  above  equation  yields  the  number  of  atoms  which  have  passed  the 
point  x  in  a  mi-infinite  slab  of  silicon  dioxide.  This  simplified  view  represents  the 
case  where  the  metal  atoms  pile  up  at  the  SiOv/Si  interface  or  diffuse  into  the 
silicon  substrate  but  do  not  interfere  with  the  diffusion  process.  The  value  for  x 
is  the  silicon  dioxide  thickness  that  we  wish  to  consider.  Figure  [3.1]  is  a  plot  of 
this  equation  versus  temperature  for  an  oxide  thickness  of  1000  A  and  a  time  of  1 
year.  It  can  be  seen  that  all  the  atoms  considered  do  not  cause  any  failures,  as  we 
have  defined  failure,  until  very  high  temperatures.  Figure  [3.2]  is  a  similar  plot, 
however  the  time  is  10  years.  We  observe  only  a  slight  increase  in  the  number  of 
neutral  atoms  that  passed  through  the  silicon  dioxide  for  the  ten  year  case  versus 
the  one  year  case.  This  result  implies  that  for  the  neutral  case,  the  solid  solubility 
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dominates  the  diffusion  process  in  the  silicon  dioxide.  The  metal  concentration 
reaches  solid  solubility  rapidly  and  no  further  diffusion  takes  place.  Fortunetly  the 
solid  solubilities  for  all  most  all  metals  is  small  enough  not  to  rule  any  out  as  possible 
candidates  for  VLSI  applications. 


2.  Ionized  Case 


Equation  ( 2.3 )  is  used  for  the  ionized  case  to  predict  a  failure.  Using  the 
same  definition  of  failure,  i.e.  10 l()  atoms /cm?  or  more  that  reach  the  interface  for 
a  given  time  and  temperature,  as  in  the  neutral  case.  We  rewrite  Eqn.  ( 2.3 )  as 

cNt  =  WDC£t  M- 

Where  Cjqt  is  the  number  of  atoms  that  reach  the  interface  per  square  centimeter. 
D  is  the  diffusion  coefficient,  Eqn.  (1.1).  C  is  the  solid  solubility,  Eqn.  (1.30).  £ 
is  the  electric  field,  106  V /cm.  T  is  the  temperature,  k  is  the  Boltzmann  constant 
and  q  is  the  charge  on  a  electron,  t  is  the  time. 

For  the  ionized  case  an  additional  criteria  for  failure  is  needed,  the  value 
of  the  electric  field.  We  set  the  field  at  a  value  of  106  V /cm.  Using  this  value  for  the 
electric  field  and  one  year  for  the  time,  Eqn.  ( 3.4 )  is  plotted  against  temperature, 
Figs.  [3.3  —  3.5],  to  reveal  the  failure  rate  for  ionized  metal  atoms  diffusing  in 
SiO 2-  To  evaluate  failures  after  ten  years  or  any  other  length  of  time,  we  simply 
multiple  the  values  in  these  figures  by  the  length  of  time  in  years. 

These  curves  depend  on  the  silicon  dioxide  thickness  only  in  the  deter¬ 
mination  of  the  effective  electric  field.  Thus  for  this  approach  to  be  valid  the  solid 
solubility  of  each  metal  must  be  reached  quickly.  The  observed  behavior  of  copper 
and  silver  support  this  assumption. 

Figure  [3.3],  the  singly  ionized  case,  predicts  the  most  failures.  Those 
atoms  that  ionize  easily  should  be  considered  very  carefully  before  using  in  VLSI 
applications.  Metais  to  the  left  of  and  including  palladium  in  Fig.  [3.3]  should  be 
safe  for  for  VLSI  applications. 
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Figure  3.4.  Doubly  Ionized  Atoms  that  have  Reached  the  Si02/Si  Interface  after 
One  Year.  Calculated  failure  rate  of  doubly  ionized  atoms  after  one 
year.  The  value  for  Na  is  -<  10r,cm~2. 
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Figure  3.5.  Triply  Ionized  Atoms  that  have  Reached  the  SiOz/Si  Interface  after 
One  Year.  Calculated  failure  rate  of  triply  ionized  atoms  after  one 
year.  Cu  and  Na  have  values  -<  105cm-2. 


The  model  presented  so  far  only  predicts  the  diffusion  and  solid  solubility  of 
atoms  in  and  through  silicon  dioxide.  There  are  two  major  considerations  that  can 
dominate  and  invalidate  the  present  approach.  The  first  consideration  is  when  the 
atom  in  question  reacts  with  the  silicon  dioxide  and  the  second  consideration  is 
whether  the  atom  reacts  with  any  other  substances  present. 

Examining  the  atoms  found  in  Appendix  A.  we  see  that  both  titanium  and 
aluminum  reduce  silicon  dioxide.  This  reaction  is  the  reason  that  aluminum  has 
worked  so  well  as  a  metallization  system.  Arsenic,  phosphorus,  boron  and  sodium 
are  know  network  formers  in  silicon  dioxide  [3.1— 3.5]  and  therefore  diffuse  by 
methods  different  than  those  presented  in  this  work. 

All  of  the  atoms  we  have  considered,  except  for  gold  and  platinum,  form 
stable  oxides.  Oxidation  of  the  diffusing  atom  will  obviously  cause  a  change  in  the 
predicted  diffusivity  and  solid  solubility.  The  data  for  copper  is  a  vivid  example. 
Another  example  is  manganese  which  oxidized  so  rapidly  that  continuity  of  the 
electric  field  was  lost  very  quickly.  Even  so  manganese  reached  levels  detectable  by 
RBS  in  the  silicon  dioxide,  in  further  support  of  our  model. 

The  reaction  of  the  diffusing  atom  with  any  substance  will  also  alter  its 
predicted  behavior.  It  is  believed  that  a  reaction  of  some  diffusing  atoms  with 
w  ier  molecules  has  a  catalytic  enhancement  on  the  diffusion  of  these  atoms.  Once 
an  understanding  of  this  reaction  is  obtained  it  can  be  easily  incorporated  into  our 
model. 

In  addition  to  the  above  considerations,  a  metal  that  is  to  be  considered  for 
VLSI  applications  must  also  have  good  adhesion,  be  easily  etched,  delineated,  form 
ohmic  contacts  with  silicon  and  be  easily  applied. 


From  this  study  in  can  be  concluded  that  Pd,  Au,  Mo,  Ta,  W,  Pt,  Ti,  and 
A1  are  potential  candidates  for  VLSI  metalization. 

It  is  possible  that  gold  should  be  eliminated  due  to  its  well  known  diffusion 
in  silicon.  Others  aspects,  such  as  adhesion,  low  contact  resistance  and  easy  chern- 
mical  processing  should  also  be  considered  before  any  final  choice  is  made. 
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APPENDIX  A 


A.  Material  Parameters 


In  this  section  the  material  parameters  used  in  the  interstitial  diffusion  model 
are  presented  as  well  as  summary  tables  of  the  copper  and  sliver  data. 


Table  A.l. 

Silver  Pile  Up  at  SOfTC. 

Diffusion  of  silver  through  silicon  dioxide  by  the  pile  up  of  silver  at  the  interface. 


A  time  time  Area  @  He?10?  (N  ). 

(sec)  (sec)  Interface  g 

(cts)  (cts)  (cm~2) 


(Nt)Ag/tlme 

(cm~2/sec) 


2700 


5400 


10800 


2700 

5400 

10800 


660 


3000 


10640 


1100  7.99x10*4 

1025  3.98x10*5 

1060  1.34x10*5 


2.96x10** 
1.19x10*2 
1.33x10*2 
1.75x10*2 
1.24x10*2 
1. 47x10*2 
1.36x10** 


21600 


24234 


1100 


2.93x10*6 


Table  A. 2 


Silver  Pile  Up  at  the  Interface  for  Various  Temperatures. 
90  minute  anneals 


temp 

°C 

Area  @ 
Interface 
(cts) 

H  SiO 

HSi  2 

(cts) 

(Nt>Ag 

(cm~2) 

275 

598 

960 

8. 8x10*4 

300 

2w/5 

1029 

3.98x1015 

335 

6167 

1081 

o  -  9.1xl014 

8.95x1015 

350 

15936 

1088 

o  =  4.3xl0l5 

1. 84xl0l6 

365 

32182 

1068 

a  «  5.65x1015 

3.06x1016 

455 

24200 

1090 

a  =  5.11x1015 

1/kT 

(eV) 


21.18 

20.25 

19.09 

18.63 

18.19 

15.94 


Table  A. 3. 

Copper  Pile  Up  at  SiO 2  Interface. 

450° C  with  a  positive  electric  field  of  5  X  IQ**  V /cm. 


4  time 
(sec) 

time 

(sec) 

Area  @ 
interface 
(cts) 

H  SiO 

HSi  2 

(cts) 

(N„) 
v  t'cu 

(cm-2 ) 

(Nt)cu/time 

(cm“2/sec) 

3600 

9765 

1194 

2.99x1016 

8.32x1012 

3600 

1.50x1016 

4.17x1012 

7200 

14525 

1139 

4.67x1016 

6.49x1012 

7200 

2.83x1015 

3.93xl0U 

14400 

15413 

1162 

4.86x1016 

3. 38x1Q12 

Table  A.4. 

Copper  Interface  Pile  Up  With  Temperature. 

Copper  silicon/silicon  dioxide  interface  pile  up  at  various  temperatures. 

temp 

°C 

Area  @ 
Interface 
(cts) 

nSi  2 

(cts) 

<Nt>c„ 

(cm“2 

1/kT  (eV) 

400 

1621 

1083 

5.48x1015 

17.24 

450 

14525 

1139 

4.67x1016 

16.05 

500 

18679 

1111 

6.16x1016 

15.01 

4.29x10^ 


350 


230 


1111 


18.63 


Table 


abl 


APPENDIX  B 


A.  General  Theory  of  Reaction  Equilibria 

For  the  case  of  constant  temperature  and  pressure  the  Gibbs  free  energy  may 
be  written  as  [B.l  ]  , 


G  =  U -TS  +  PV  (B.l) 

since  it  is  at  a  minimum  with  respect  to  the  proportions  of  the  reactants.  The 
incremental  Gibbs  energy  may  be  written  as  (B.l], 

dG  =  Y,  HjdNi  =  dU  (B.2) 

3 

where  fij  is  the  chemical  potential  of  species  j  and  dNj  is  the  change  in  the  number 
of  molecules  of  species  j. 

The  change  in  the  Gibbs  free  energy  in  a  reaction  is  seen  to  be  closely  related 
to  the  chemical  potentials  of  the  reactants  [B.l].  In  equilibrium  this  change  must 
be  zero  for  constant  P  and  T. 

dG  =  ^2  HjdNj  =  0  (B.3) 

3 
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The  change,  dNj,  in  the  number  of  molecules  of  species  j  is  proportional  to 
the  number  of  moles  of  species  j  as  determined  by  the  balanced  chemical  equation 
of  reaction.  Let  Uj  denote  this  coefficient  and  let  Aj  denote  the  chemical  compound 
of  species  j.  Then  the  balanced  chemical  equation  can  be  written  as 

E  =  0  (B.4) 

3 

We  can  write  dNj  as: 

dNj  =  UjdN ’  (B.5) 

where  dN ’  is  the  increment  in  the  number  of  times  the  reaction  takes  place  [B.l]. 
Combining  Eqn.  (B.3)  and  (B.5)  results  in: 

«  =  £««'=  0  (B-6) 

3 

simplifying 

E  W  =  0  (B.7) 

3 

Equation  ( B.7 )  must  be  satisfied  to  be  in  equilibrium  at  constant  temperature 
and  pressure. 

B.  Law  of  Mass  Action  and  Equilibrium  for  Idea  Gases 

If  we  assume  that  each  constituent  behaves  like  an  idea  gas,  the  chemical 
potential,  fij,  of  identical  polyatomic  molecules  is  found  as  follows  (B.i  ]: 


Each  molecule  has  rotational  and  vibrational  motions  so  that  the  total  energy, 
Et,  of  a  molecules’  orbital  t  is  the  sum  of  two  parts, 


Et  —  Etnt  +  En 


where  refers  to  the  internal  degrees  of  freedom  and  En  to  the  translational 
motion  of  the  center  of  mass  of  the  molecule.  For  translational  motion  it  can  be 
shown  [B.l]  that 


fi2  /7rn\2 


Jn  2m  V  L  ) 


where  n  is  the  quantum  number  of  the  translational  orbital  and  L  is  the  is  integer 
times  half  a  wave  length.  E{nt  includes  the  vibrational  energy  of  the  molecule  on 
the  force  constant  and  the  rotational  energy  of  the  molecule  [B.l]. 

In  the  classical  regime  the  probability  that  a  given  translational  orbital,  n,  be 
occupied  is  always  very  small  in  comparison  with  one.  When  we  write  the  grand 
sum  of  the  orbital  in  this  regime  we  may  neglect  terms  in  the  X2  and  higher  powers 
of  X,  because  such  terms  correspond  to  occupancy  of  the  orbital  by  more  than  one 
molecule. 

Accordingly,  the  grand  sum  for  the  system  of  all  orbitals,  t,  for  which  the 
translational  quantum  number  is  precisely  n  and  for  which  the  internal  quantum 
number  i  assumes  all  possible  values  is 


z  —  1  +  X  ^2  e  (E">t+E’')/kT 


(H.  10) 


since  we’re  assuming  the  classical  regime  we  may  factor  e  E''/kT  to  obtain 
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z  =  1  +  Xe 


-En/kT  ' 


-Et/kT 


(B.11) 


Note  that  X  is  the  absolute  activity,  X  =  .  We  define 


^  =  Eint,kT  (B.12) 

t 

as  the  partition  function  of  all  the  internal  degrees  of  freedom. 

So  equation  (B.ll)  becomes 


z  =  l  +  \zinteEn!kT  (B.13) 

The  probability  that  the  translational  orbital,  n,  is  occupied  irrespective  of  the 
state  of  internal  motion  of  the  molecule,  is  given  by  the  ratio  of  the  term  in  X  to 
the  grand  sum  z. 


f{En)  = 


■\zinte-E»/kT 


1  +  Xz,nte-En/fcT 
for  the  classical  regime  which  holds  when  f{En )  «  1 


^inte~En,kT 


(B.14) 


If  we  now  compare  this  result  to  that  of  an  idea  monatomic  gas  [B.l  ]  we  note 
that  the  only  difference  is  the  addition  of  the  z,-nj  term  and  X  becomes: 


X  = 


(B.1S) 


and  thus 


* 


i 


H 

« 

t 

V 


t 


2 


N  / 

'27Th2\5  ,  % 

Mj  =  kT 

A 

^mjkTj  ln(Z%nt) 

(B.16) 


If  we  let  Cj  =  N  jV  denote  the  concentration  of  molecules  of  species  j  we  can 
rewrite  equation  (B.16)  as 


Hj  =  kTlnCj  +  -kTln 

u 


f  2ir  H2  \ 
\mjkT  ) 


+  E, 


int 


(B.17) 


Note  that  2?inj  includes  all  internal  degrees  of  freedom  regarding  vibration, 
rotation,  electronic  excitations  and  all  the  nuclear  orientations. 


This  result  for  the  chemical  potential  is  the  sum  of  the  term  in  the  logarithm 
of  the  concentration  and  a  term  that  is  a  function  only  of  temperature: 


Mi  =  kTlnCj  +  kTXj(T) 


(B.18) 


where  we  define 


Xj 


= -4^) 

2  \rrijkT  J 


+ 


E, 


int 


kT 


(B.19) 


Note  that  the  internal  free  energy  E{nt  is  an  additive  term  in  the  chemical 
potential  of  the  jth  chemical  component. 


The  equilibrium  condition  ^iHj  —  0  now  becomes 


2Z  viVj  =  kT  J2(U'lnCJ  +  U'Xj)  =  0 


(B.20) 


or 
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APPENDIX  C 


A.  Diffusion 

Diffusion  is  the  random  molecular  transport  of  material  to  relieve  a  non¬ 
equilibrium  condition.  This  transport  is  driven  by  a  concentration  imbalance  where 
on  the  average  the  areas  with  a  higher  concentration  of  material  loses  or  trans¬ 
ports  material  to  areas  of  lower  concentrations.  Additional  driving  forces  such 
as  electromagnet  fields  or  perhaps  gravitational  fields  can  aid  or  hinder  diffusion 
depending  on  their  direction. 

In  this  section  the  diffusion  equations  for  positively  and  negatively  charged 
particles  will  be  derived.  First  we  will  derive  the  flux  for  a  charged  particle  and 
the  transport  equation  which  determines  their  distribution  in  a  solid.  Then  we  will 
present  solutions  to  the  transport  or. diffusion  equations  under  various  conditions 
and  assumptions. 


1.  Flux 

Figure  [  1  ]  represents  the  potential  hills  formed  by  the  atoms  of  a  crystal. 
These  potential  hills  impede  the  motion  of  any  charged  particles  within  the  crystal. 
The  height  of  the  potential  barrier  W  is  typically  of  the  order  of  electron  volts 
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CL 


x-a  x  x+a 

DISTANCE 

Figure  C.l.  Potential  Distribution'.  Model  of  ionic  motion  within  a  crystal  with 
zero  electric  field. 

[C.l  — C.3]  in  most  materials.  The  distance  between  successive  potential  barriers, 
a,  is  of  the  order  of  a  lattice  spacing  which  is  several  angstroms. 

If  a  constant  electric  field  is  applied,  the  potential  distribution  as  a  func¬ 
tion  of  distance  will  be  tilted  as  shown  in  figure  [2].  This  will  help  positive  ions 
move  to  the  right  and  hinder  their  motion  to  the  left.  With  the  field  applied  the 
flux,  J,  ( J  =  number  of  any  species  passing  through  a  unit  area  in  a  unit  time) 
at  x  will  be  the  average  of  the  fluxes  at  position  ( x  —  |)  and  at  (z  +  |).  Using  the 
notation  shown  in  figure  [2], 


r,_,  _  (/I  -  J 2)  +  (/ 3  -  / 4) 
j[x)  —  - 


(C.l) 


Consider  the  component  J.  It  will  be  given  by  the  product  of  (l)  the 
density  per  unit  area  (normal  to  the  direction  of  the  flux)  of  particles  at  the  potential 
valley  at  (z  —  a),  (2)  the  probability  of  jumps  of  any  of  these  particles  to  the  valley 
at  z  and  (3)  the  frequency  of  attempted  jumps,  i/. 


Jl  =  aC{x-a)e'3’lw-*a£'>is 
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(C.2)  | 

i 


V  •«  .■.vv.f  jy. -'J 


Figure  C.2.  Potential  Distribution.  Model  of  ionic  motion  within  a  crystal  with  an 
applied  electric  field. 

where  aC(x  —  a)  is  the  density  per  unit  area  of  particles  situated  in  the  valley  at 
(i  —  a)  [C.2  ,  C.4]and  the  exponential  factor  is  the  probability  of  a  successful  jump 
from  the  valley  at  ( x  —  a)  to  the  valley  at  x  [C.4,  C.5  ].  Note  the  lowering  of  the 
hairier  due  to  the  electric  field,  £ . 

Likewise 


J2  =  aC 7(x)e_^u'+5of>i/ 


(C.3) 


/3  =  aC{x)e-Mw-2a£K  (C.4) 

JA  =  aC{  x  +  a)e-ib(w+*a£)  u  (C.5) 


Substituting  these  expressions  for  J 1,  J 2,  J 3,  and  J 4  into  equation  ( C.l ) 


making  the  assumption  that  [C.6], 


which  is  the  transport  equation  in  one  dimension.  If  we  expand  our  thinking  to 
three  dimensions  the  transport  equation  becomes 


dC{x) 


=  —div  J[x ) 


(C.16) 


Combining  the  transport  equations  and  our  equation  for  flux,  equation 
(C.10)  results  in 


_  -±\-Dcosh 
dt  dx 


(qoS\ 

\2kTJ 


dC{x)  2D  .  ,(qa£ 

. ■~—L  +  — stn/i  - 

dx  a  V  2 


If  we  now  assume  that  D  and  c  are  independent  of  x  and  £  «  H  * 
may  write  equation  (C.17 ),  the  one  dimensional  diffusion  or  transport  equation,  a 


3C(x)  ^  dC2(x)  dC(x) 
dt  dx 2  ^  dx 


Likewise  for  negative  species  under  the  same  assumptions 


(C.18 


dC{x)  _  dC\x)  dC{x) 
dt  dx 2  P  dx 


(C.19 


B.  Solutions 

Before  solving  the  one  dimensional  diffusion  equation  let  us  reiterate  the  as 
sumptions  imposed  so  far  and  comment  on  their  validity: 
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(1)  Diffusion  only  in  one  dimension,  normal  to  the  surface. 

(2)  The  electric  field  and  the  diffusion  coefficient  are  independent  of  time  and 
direction. 

(3)  e «  a?. 

(4)  Solid  is  homogeneous. 

Assumption  one  is  reasonable  since  we  are  considering  the  diffusion  of  charged 
particles  and  thus  diffusion  is  enhanced  parallel  to  the  direction  of  the  electric  field. 
The  electric  field  is  independent  of  time  and  direction  if  the  solid  is  locally  uniform. 
The  assumption  that  D  is  independent  of  time  and  direction  requires  a,  /x  and  w  be 
independent  of  time  and  direction.  This  is  true  if  the  solid  is  locally  homogeneous. 

To  see  if  £  <C  we  assume  a  =  i-0  A  which  is  a  conservatively  high  value 
for  a  lattice  spacing.  Then 


1723.40T  V/cm 


(C.20) 


where  k  is  the  Boltzmann  constant  and  q  is  the  charge  of  an  electron. 

Table  1  shows  for  different  values  of  temperature.  We  see  that  fields  of 
104  —  105  V/cm  are  probably  small  enough  to  make  assumption  three  valid. 

Table  C.l. 

Estimates  of  i  <£. 


Temp  °K 


5.170x105  v/cm 
1.034xl0&  V/cm 
1.206xl06  V/cm 


Assumption  four  is  valid  since  our  samples  are  Czochralski  grow  silicon  with 
thermally  grown  silicon  dioxide  and  is  therefore  homogeneous,  at  least  locally. 

Using  these  assumption  we  solve  the  diffusion  equation  in  one  dimension  for 
positive  particles,  equation  (C.18),  for  the  following  boundary  conditions. 

(1)  C{0,t)  =  Ca. 

(2)  (7(oo,  t)  =  0. 

(3)  (7(1,0)  =  0. 

where  Ca  is  the  solid  solubility  of  the  particle  in  the  silicon  dioxide.  The  solution 
can  be  shown  to  be  [  C.7 — C.9  ]. 


C(x)  = 


q£x 

+  e'ET"  erfc 


(C.21) 


Figure  [4]  is  a  plot  of  this  equation  for  300°C,  1.5  hours,  £  ==  5.2  X  104V/cm 
and  a  solid  solubility  of,  Cs,  of  3  X  1018cm~3.  Different  diffusivities  are  shown  in  the 
figure  to  allow  a  comparison  with  the  silver  data  presented  eatlier.  This  solution  to 
the  diffusion  equation  models  reality  when  the  substance  the  particles  are  diffusing 
through  is  semi-infinite.  As  a  further  example  figure  [  5  ]  is  a  plot  of  equation  ( C.21 ) 
when  T  —  100°C,  t  =  one  year,  £  =  106  V/cm,  Ca  =  3  X  1018cm-3,  and  the 
diffusivity  takes  on  the  values  shown  in  the  figure. 

As  an  approximation  to  the  case  when  the  diffusing  particles  pile  up  at  the 
SiOi/Si  interface,  equation  (C.18)  may  be  solved  with  following  boundary  condi¬ 
tions. 


(1)  C(0,  *)  =  (?,. 

(2)  C(x,  0)  =  0. 


where  a  is  the  film  thickness. 


The  solution  to  these  conditions  is  [C.7  —  C.9], 


C{x,t)  = 


e  e2ki  stn 


(C.22) 


where  the  /?n’s  are  the  positive  roots  of 


tan 


(C.23) 


Figure  [6]  is  a  plot  of  this  equation  for  300° C,  1.5  hours,  £  =  5.2  X  104V /cm, 
C9  —  3  X  1018cm~3,  and  a  films  thickness  of  5700 A.  Figure  [7  J  is  a  plot  of  equation 
(C.22)  with  T=  lOO0^,  £  =  10 *V/cm,  C9  =  3  X  1018cm~3,  a  —  5700  A,  and  a 
time  of  one  year. 


022*10  crrr/sec 


D=  10  cm^/sec 


D«  10‘15cm2/sec 


■0  =10*16cm2/sec 


TEMP:  300°C 
TIME  :  1.5  hrs. 
Cs  =  3  *1018cm*3 
£  *5.2  *  10*  V/cm 


DIST  (k  A) 


Figure  C.4.  Solution  to  the  Diffusion  Equation  for  a  Semi-  infinite  Slab.  300°C', 
1.5  hours,  if  =  5.2  X  10 4  V/cm  and  Ca  =  3  X  1018cm-3. 


D>10  wcm  /sec 
0  *  10’2®cm2/sec 
D=  5  x  I0'21cm2/sec 
D*  10'21cm2/sec 


TEMP:  100°C 
TIME  :  1  YR. 
Cs  =  3xio1acm 
£= 10®  V/cm 
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Figure  C.5.  Solution  to  the  Diffusion  Equation  for  a  Semi-  infinite  Slab.  100o(7, 
t  =  one  year  ,  C  =  10 G\//cm  and  CB  =  3  X  10,8cm-3. 


0  1  2  3  4  5  6 
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C.6.  Solution  to  the  Diffusion  Equation  for  a  given  film  thickness.  300° C, 
1.5  hours,  i  =  5.2  X  10 4V/crn,  C3  =  3  X  1018cm-3  and  a  film 
thickness  of  5700  A. 


DIST  (k  A) 

C.7.  Solution  to  the  Diffusion  Equation  for  a  given  film  thickness.  100°C’, 
t  —  one  year  ,  £  =  106 V/cm,  Ca  —  3  X  1018cm-3  and  a  film 
thickness  of  5700  A. 
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